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Within a minimal model for the iron-based superconductors in which itinerant electrons interact 
with a band of local moments, we derive a criterion for sign-reversing superconductivity. We show 
generally that even if the inter-band and intra-band interactions are repulsive, sign-reversing super- 
conductivity obtains as long as the inter-band scattering dominates. This effect can be modeled as 
arising from an internal Josephson link between the Fermi surface sheets. Although the experimen- 
tally observed isotope effect is large, a = 0.4, we show that it is consistent with a non-phononic 
mechanism in which it is the the isotope affects results in a change in the lattice constant and as a 
consequence the zero-point motion of the Fe atoms. 
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Although the history of superconductivity in iron- 
based materials is shorljii^i^, a consensus seems to have 
been reached that phonons are not the efficient cause 
of the pairing mechanism. For example, density func- 
tional calculations'* predict that for LaFeAsOi^ajFa;, the 
dimensionless electron-phonon coupling constant is A = 
0.21 and the logarithmic-averaged frequency is 206K. 
In conjunction with the standard BCS expression for 
Tc, these values lead to a transition temperature of 
0.8K which is too low to explain 26K superconductiv- 
ity. Experimentally, the bare phonon density of states^ in 
LaFeAs0i_2;Fa; (hereafter 1111) is in perfect agreement 
with the density functional results^ and hence corrobo- 
rates that phonons are too weak to explain superconduc- 
tivity. As a result, mechanisms, such as spin fluctuations, 
which exploit the proximity of the superconducting phase 
to the antiferromagnet with (tt, 0) order— in the parent 
material have risen to the fore^. 

Consequently, it has come as a surprise that the iron- 
pnictdes have a sizeable isotope effeclj^. When ^^Fe is 
replaced by '^''Fe in Bai_2:K2,Fe2As2 (hereafter 122), the 
transition temperature for the magnetic order at x = 
and superconductivity at x = 0.4 is changed by 1.39% 
and 1.34%, respectively. Expressing the mass depen- 
dence of Tc using a = — din Tc/d In M, leads to a value 
of a = 0.4 for the observed change in Tc. A full iso- 
tope effect in standard low Tc materials corresponds to 
a = 0.5^. Hence, the pnictides at the doping levels stud- 
ied have an isotope effect indicative of elemental super- 
conductors. The origin of this effect is not known but 
certainly presents a challenge for purely electronic mech- 
anisms of superconductivity in which a is supposedly 
negligible. 

In this paper, we address the apparent contradic- 
tion between the irrelevance of phonons to Cooper pair- 
formation and the sizeable isotope effect in the pnic- 
tides. Aside from the isotope effect, any pairing mech- 
anism must resolve the high critical temperature and 
the competition and between magnetism, that is, the 
(tt, 0) antiferromagnet, and superconductivity. Within 
a two-band modeU^ for the Fe-As layer in which local 
moments in one band interact with itinerant electrons 



in the other, we show that interband scattering medi- 
ated by magnons leads to a high critical temperature. 
Our work provides a microscopic framework for apply- 
ing the well-known result^' that sign-reversing pairing 
does not require attractive inter nor intra-band interac- 
tions in two-band models, as in the case of the pnictides. 
Second, we show that although the pairing in this model 
is electronic in nature, an isotope effect exists as a result 
of the sensitivity to the lattice constant. An a = 0.4 re- 
quires a zero-point motion of the Fe ion to be 6.2meT^, 
in agreement with experimentally known^^ values. 

An important fact about the pnictides is that they 
exhibit a (tt, 0) antiferromagnet but they are metallic, 
nonetheless. Metallicitiy in the presence of such order 
implies that some states at the chemical potential re- 
main ungapped. Whether the magnetism arises from lo- 
cal or weak-coupling physics is currently not resolved^-. 
However, models in both of these extremes are prob- 
lematic. For example, pure spin models such as the 
J1-J2 modeliiii^ii^ which rely on a fine-tuning of the 
nearest and next-nearest neighbor exchange interactions 
are clearly incomplete as they describe insulators. At 
the other extreme, weak-coupling scenarios based on 
nestingiiii^ii^ are also problematic because the mono- 
clinic distortion"^' 1^ that precedes the magnetically or- 
dered phase in the pnictides shifts each half-filled band 
away from the perfect nesting condition by an amount 
related to the crystal field splitting. Regardless of which 
model is used, it must at least account for the experimen- 
tal fact that local magnetic correlations exist above the 
ordering temperature as evidenced^S by the temperature 
dependence of the peak-to-peak linewidth and g factor in 
electron spin resonance studies on LaFeAsOi-ajFa,. The 
simplest scenario that is consistent with this physics is 
a two-orbital model in which an itinerant band of elec- 
trons hybridizes with a set of local moments residing on 
the other level. In a previous paper**^, we have shown 
that the multi-orbital physics of the Fe-As interactions 
reduces to an effective two-orbital model. The itinerant- 
localized dichotomy arises from the difference in the p-d 
character of the two levels. While the p-d hybridiza- 
tion in this model appears as a tuning parameter, it was 
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foundlS that a hybridization of 0.8eV was needed to ex- 
plain the magnitude of the magnetic moment. This value 
of the hybridization is consistent with the only available^i 
experimental estimate. A further key prediction of this 
work is that aside from the moment lying in the a — & 
plane, a residual moment lies along the z— axis with a 
magnitude of 0.06/^5, as found experimentally^^. 

In light of the success of this model in describing the 
parent magnetically ordered state, we adopt it here to 
investigate possible superconducting instabilities. The 
model can be written as a spin-fermion Hamiltonian^^i^l, 

H = H, + H, + Hsi (1) 

He = £k,g^c]^.^Ck,(T 

= Ji ^ Si • Sj + J2 51 Si • Sj 

n.n 71. n.n 

Hsf = J ^ cl^CTapCi^iS ■ Si. (2) 




FIG. 1: a) Fermi surfaces and Cooper pairs (CP) of Iron- 
based SC on the folded Brounian Zone. Two hole pockets 
(a,/?) on r point and two electron pockets (71,72) at M point 
are determined by dispersion relation — ^ and -B^^q = ^ 
respectively. Of two possible interband CP scattering process, 
(k I, — k [) scatted to (k I, — k and (k — k J,) scatted 
to (k' -I- Q J,, — k' -|- Q 1)5 only the first case are momentum 
conserved, b) Interaction vertex mediated by magnon. c) 
Attractive interaction mediated by interband scattering. 



in which describes the itinerant electrons, the lo- 
calized electrons which yield the magnetism and i/gf the 
spin-fermion interaction between the two sets of elec- 
trons. Here ek is the band structure which yields the 
Fermi surface of the non-interacting system, cj^. ^ creates 
an electron with momentum k and spin u, and Si repre- 
sents the spin on site i. We have retained only the nearest 
and next-nearest spin exchange interactions, Ji and J2, 
respectively, as experiments and theor y^^'^^i^^ indicate 
that interactions beyond these are negligible. Although 
Ji and J2 will be treated as phenomenological param- 
eters, it important to note that they are both directly 
proportional to the p-d hybridization. That is, if the p-d 
hybridization is ignored, as advocated recently2£, both 
vanish. 

We first perform a Holstein-Primakoff transformation 
on the spin-part of the Hamiltonian to obtain an effective 
description of the spin excitations on the (tt, 0) antifer- 
romagnet. To order l/S*, we obtain^!, 

Hs = Ci + S'^[ylkaj^ak + ^(Skfl^a^^ + B*kaka-k)] 
Ak — 4 J2 -f- 2 Ji cos kx 
-Bk — 4 Ji cos ky + 8J2 cos kx cos ky 



(3) 



When diagonalized using 6k — cosh^kik — sinh^koLk ' 
this Hamiltonian yields a dispersion for the magnons of 
the form 



= C2 + 5^c^k6|,&k 



c^k = S^/Al-Bl 



(4) 



where Ci{i = 1,2) are all constant. It is now convenient 
to express the Fourier transform of our original Hamilto- 
nian in terms of these effective spin excitations. By us- 

y^j:^e-'^'a^, and 



ing the Fourier transform cig 
Cma = yWEke^^^^+'^/^ka whcrc a 



and m index the up and down spins respectively, we can 
rewrite the Hamiltonian^!, 



H 



sr 



(5) 
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-J^/2S ^(5(k'-k + q 

k,k',q,P v=± 

/,.,,(k,k',p)4r;<x«q+'^-c- 

i?ff = JV2S [^o-cj^cio^ {a\a\ - (a\a\ 



P) 



(!■ 
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in terms of the the spin normal modes, Ok and the elec- 
tronic degrees of freedom, d^^, which are linear combi- 
nations of ak(T and /Jka- The parameters in this effective 
Hamiltonian are defined by 



AS - (ek - ek+q) /2 , ^k ^ 
Et = (ek + ek+Q)/2±i?k. 



(6) 
(7) 



(a, 0) and 1 



where P is the crystal reciprocal lattice, Q = (0,7r/a) 
in the unfolded Brillouin zone represents the wave- vector 
for the magnetic order, /i^^i,/ (k, k', P)^ is a function of 
P, 6'k and 6'k' where sin6'k = ~ J/E^) /2. The full 
details, which can be found elsewhere^!, of /i,^^ are not 
particularly relevant as only its magnitude appears in the 
final equations. We emphasize that ek is the bare spec- 
trum of the electron in the itinerant level. However, as 
a result of the antiferromagnetic correlations, the spec- 
trum of the itinerant electrons is split into two branches 
E^{k) which satisfy £'='=(k) = E^{'k + Q) in the ordered 
state. However, in the doped case, in which it is only 
the remnant spin fluctuations which persist, this perfect 
nesting condition is destroyed. 

There is of course an extreme similarity between the 
various contributions to Eq. ^ and the electron-phonon 
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Hamiltonian used in BCS. In this case, the phonons are 
replaced by magnons and the interactions in involve 
spin flips. It is from this term that the dominant interac- 
tions arise. Treating this term perturbatively to second 
order, we obtain the amplitude 



(2j25V(q)|/,,,.(k,k',P)|^ 



c.2(q) 



(8) 



for the scattering of an electron pair with momenta 
(k, — k) to one with momenta (k', — k'). In Eq. ([8]), €e{^) 
is the dielectric constan t^^'^^ and q = k' — k. Regardless 
of the momentum transfer, the interaction arising from 
magnon scattering is always repulsive as a result of the 
minus sign in front of Eq. ([8]) . This crucial difference^^ 
with the phonon-mediated interaction in BCS arises from 
the spin-flip nature of the magnon scattering. 

Nonetheless, sign-reversal superconductivity is still 
possible even though the scattering processes are 
repulsive. To establish this we consider the sep- 
arate amplitudes for inter-band (between the M 
and r points) and intra band (within the M or F 
points) scattering processes. Intraband scattering 
dominates in the long wavelength limit. Rewriting 
qx — qcosa and qy = gsina, we have that ^(q) = 

2Sq^{2J2 + Ji) [(2 J2 - Ji) cos2 a + (2J2 + Ji) sin^ a] = 

2SqF(a). In this limit, qte(c{) sa q A- qs where 
qs = 2ne^N{EF) is the inverse Fermi-Thomas screening 
length and N{E-p) is the density of states around 
the Fermi surface. For the intraband scattering, 
^k' ~ -^k ~ 0. Consequently, the magnon- mediated 
interaction is of the form, 



^i(k,k') = 



(J2)l^..(k,k^p)|- 

qsF{a) 



(9) 



where ii denotes ee for intra-band scattering within an 
electron pocket and hh for scattering within a hole 
pocket. We compare this amplitude to that for inter- 
band scattering. The only interband scattering process 
that is momentum conserving involves the scattering of 
an electron pair from the hole pocket at F to two different 
electron pockets such that the final pair has momentum 
(k' -I- Q i, — k' — Q t). For the interband scattering, we 
still have E^, — E^ « but because the momentum car- 
ried by the magnon is Q -f q, quantitative differences. 



(J2)l/....(k,k^p)|- 

e(Q + q)|q|G(a) 



(10) 



arise with the intra-band scattering amplitude. Note the 
presence of the inverse \q\ dependence. Here, G{a) = 
\/2Ji(2J'2, -|- Ji). Since Q = (7r,7r), it is no longer ap- 
propriate to take the long- wavelength limit of the dielec- 
tric function. As a result of the q~^ divergence of Eq. 
(|10p . the inter-band interaction is much larger than the 
intra-band interaction for small momentum transfer |q|. 
This is the first demonstration from the microscopies that 
inter-band scattering dominates the pnictides. 



In the limit that the inter-band scattering dominates, 
the BCS gap equations admit sign-reversing pairing en- 
tirely from repulsive interactions. This can be seen by 
solving the gap equations for the 



A 



E 



Vhh Kill 



(k,k') 



FS 



for a 2-band system. A crucial difference with earlier 
work is the minus sign on the left-hand side of Eq. (|lip . 
If we assume all Fermi surfaces are circular, then inte- 
grating over the Fermi surface is equivalent to multi- 
plying by the DOS. Defining the dimensionless interac- 
tion strength Ac = Ycc^ciEY-), Ah = VhhA'^h(-BF) and 
Aoh = T4h\/A'c(-£'F)-/Vh(-£'F), we obtain the maximum 
positive eigenvalue. 



A 



(Al, + Ae)/2+A/(Ah-Ae)V4+A2i^, (12) 



and the corresponding eigenstate (Aj^,A^) oc 

(Ach, — (A + Ah))'^. We see clearly that for A > 0, the or- 
der parameter on the hole and electron Fermi surfaces has 
different signs only if Aoh is positive and exceeds all other 
couplings. This is our principal conclusion and shows 
that the multi-band structure of the pnictides leads nat- 
urally to s± pairing symmetry. The critical temperature 
is given by 



Te = 1.14afexp(-l/A), 



(13) 



where ilgf is the counterpart of the Debye frequency for 
magnons, which is estimated to be 1500K for J\ = J2 — 
SOOK^^.. Coupled with a large A as a result of the domi- 
nance of Aoh, the transition temperature can be sizeable. 

A simple physical argument is instructive here to de- 
lineate why pairing can obtain from repulsive inter-band 
interactions. Consider the second-order process shown in 
Fig. llJ:). A Cooper pair at a hole pocket is scattered 
into a Cooper pair at an electron pocket and then scatted 
back to the hole pocket. Regardless of the sign of the in- 
terband scattering, the resultant matrix element for this 
process 



[/(k,k') 



/da;/d3pyoh(k,p)T/oh(p,kO 
(cj + i(5 - Cp+q) (-r^ - - Z(5 - e-p-q) 



,3 yoh(k,p)Kh(p,k') 



d'^p 
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< 0, 



Cp+q 



p Q 



(14) 



is always attractiv o^i— . Here, — = ^(^k + c-k) is the 
energy of two electrons on the hole pockets. The key 
point here is that as long as electron and hole pockets 
are separated in energy, an attraction develops through 
an exchange of two electrons rather than single-electron 
hopping. In this sense, the interband scattering can be 
viewed as an internal Josephson junction with coupling 
Jtun cos((?!)i — ^-2) where 01 and ^2 are the phases of the 
order parameter on the two Fermi surfaces. 
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A striking feature of the pnictides is the relative in- 
sensitivity of to doping. Within our model, this is 
naturally explained because the key factors that deter- 
mine Tc are the magnon energy and the coupling J. Since 
J is simply determined by the localized carriers, doping 
has a relatively negligible effect. Predominantly, doping 
changes the density of itinerant electrons as is confirmed 
experimentally'^^ . 

What about the isotope effect? In the phonon- 
mediated case, the isotope effect appears through the de- 
pendence on the Debye frequency. Because the courter- 
part of the Debye frequency, Qgf, for spin- fluctuation 
pairing is not related to the ion mass, there should be 
no such effect. However, beside changing the Debye 
frequency, isotopic substitution also modifies the lat- 
tice constant a^. Smaller values of a will increase the 
electron-electron interaction, thereby enhancing the spin 
flucuations and, as a consequence, a higher Tc. Exper- 
imentally, the critical temperature of the 122 materials 
is sensitive to pressure^ and hence consistent with this 
explanation. Experimentally, and TgDw decrease at a 
rate of AT^/AP « 0.22K/kbar for (Bao.55Ko.45)Fe2As2 
and ATsDw/AP w l.OK/kbar for BaFe2As2^^. Conse- 
quently, it is reasonable to consider a harmonic toy model 



in which the isotopic effect arises from a modulation in 
the lattice constant. Within such a heuristic model, the 
isotopic parameter is given by 



ATjTc ^ 1 
' AM/M 2 



sinh(a;) 



> 0, 



(15) 



where x = Ticu/kBT. For Bao.6Ko.4Fe2As2, Tc ~ 36K 
and a — 0.4^. Explaining this shift requires a zero-point 
energy of iron of roughly 6.2meV. For the 122 materi- 
als, the zero-point energy of iron is not known. However, 
for iron metal, nuclear resonant inelastic x-ray scattering 
place the zero-point energy at BmeV'^^. Hence, our esti- 
mate used here to obtain the isotope shift is reasonable. 
This kind of estimate is also applicable to the cuprates in 
which a positive isotope effect is seen experimentally^. 
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